Dynamics of a bouncing dimer 
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We investigate the dynamics of a dimer bouncing on a vertically oscillated plate. The dimer, 
composed of two spheres rigidly connected by a light rod, exhibits several modes depending on 
initial and driving conditions. The first excited mode has a novel horizontal drift in which one end 
of the dimer stays on the plate during most of the cycle, while the other end bounces in phase with 
the plate. The speed and direction of the drift depend on the aspect ratio of the dimer. We employ 
event-driven simulations based on a detailed treatment of frictional interactions between the dimer 
and the plate in order to elucidate the nature of the transport mechanism in the drift mode. 

PACS numbers: 05.45.-a, 45.70.-n, 47.52. +j 
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The shape of grains has a fundamental impact on static 
and dynamic properties of granular materials. Packing of 
granular material is affected by the relative orientation of 
grains 0, the statistics of avalanching is apparently dif- 
ferent in regular sandpiles and piles of elongated rice Q , 
and strongly shaken rods exhibit self-organization which 
leads to formation of large vortices 0, y[ . Thus the phe- 
nomenology of asymmetric grains is very different from 
that of spherical beads and irregularly shaped sand. 

In Ref. 0, using a combination of theory and exper- 
iments in a quasi-2D geometry (vibrated rods in an an- 
nulus), we showed that the motion of rods is caused by 
the non-eccentric frictional collisions between the rods 
and the vibrated plate. These collisions are translated 
into collective motion by multiple non-elastic collisions 
among the rods during the flights between rod-plate col- 
lisions. In this paper, we focus on an even simpler case of 
a single anisotropic "grain" on an oscillated plate in or- 
der to study mechanisms of transformation of vibrational 
excitation into horizontal transport in the absence of ex- 
ternally imposed asymmetry. Recent studies of thermal 
ratchets have addressed the more straightforward situa- 
tion where there is an imposed left/right asymmetry in 
the substrate 0- 

In contrast with our previous work with thin rods, here 
we study the dynamics of a rigid dimer consisting of two 
heavy spheres connected by a light rod. Unlike the rods 
case, the points of contact with the plate are easily deter- 
mined even when the dimer is horizontal at the time of 
collision (in this case, there are two well defined contact 
points). We also note that the oscillated dimer can serve 
as a natural generalization to the much studied system 
of a particle on a vibrated plate which is a paradigm for 
period doubling and transition to chaos pfl. As we show 
here, in addition to those regimes, the oscillated dimer 
displays spontaneous symmetry breaking and ratchet-like 
transport in a system without externally imposed asym- 
metry. 

Figurc^shows a photo and a schematic diagram of the 
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FIG. 1: (a) Photo of a dimer bouncing on the vibrated plate, 
(b) A sketch of the system. 



system. The experiments reported here are performed 
with dimers fabricated with steel spheres with a diam- 
eter of 9.5 mm, and a glass rod with a 4 mm diameter. 
The length of the rod is varied to change the aspect ra- 
tio A r (defined as the ratio of the length of the dimer 
along the symmetry axis and the diameter of the ball) 
between 2 and 10. The dimers are placed on a plate be- 
tween two vertical glass walls which are rigidly attached 
to the plate and separated by 10 mm in order to confine 
the motion to a vertical plane. The container is vibrated 
with an electromagnetic shaker and a sinusoidal wave- 
form generator. The oscillation frequency / and magni- 
tude r = 4ir 2 Af 2 /g, where A is amplitude of plate oscil- 
lation, and g is gravitational acceleration, can be varied 
independently. The period of plate oscillation T = l/f. 
The motion of the dimer and the plate is recorded using a 
Phantom v5.1 camera at 2000 frames per second using a 
standard reflection lighting configuration and a centroid 
algorithm which yields positions to within 20 /im. 

For r > 1, complicated irregular hopping and tumbling 
motion is observed for most initial conditions. Here, we 
focus on the weak excitation regime T < 1. Clearly, for 
r < 1 the dimer may remain on the plate. However, if 
the dimer is dropped onto the plate, three main modes of 
persistent motion can be observed depending on initial 
conditions. The vertical coordinates of the tips of the 
dimer in the laboratory frame of reference are plotted 
in Fig. In the highest energy mode, the tips hit the 
plate out of phase at every other period of oscillation 
[see Fig.|2Ia)]. We call this motion the flutter (F) mode. 
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FIG. 2: The vertical position of the left (o) and right (A) 
tips @ of the dimer as a function of time for the lowest energy 
modes observed in the experiment, (a) The flutter (F) mode, 
(b) The jump (J) mode, (c) The drift (D) mode. The solid 
curve corresponds to the vertical position of the plate (A r — 
3.9, / = 50 Hz, r = 0.9.) 



In the next lower energy mode, both tips collide with 
the plate once every cycle [see Fig. Efb)]. We call this 
the jump (J) mode. In this mode the dimer does not 
have significant angular momentum and resembles the 
first excited mode of a single particle bouncing on an 
oscillated plate. The lowest energy mode is shown in 
Fig. EKc). As can be noted from the figure, one of the 
ends of the dimer hits the plate once every cycle while 
the other appears to stay on the plate throughout the 
cycle. In addition to the non-trivial vertical motion, a 
persistent horizontal drift is observed and therefore we 
call it the drift (D) mode ||. Although the excitation 
of a particular mode depends sensitively on the initial 
conditions, the drift mode can be excited easily by lifting 
one end by a fraction of the sphere diameter and releasing 
it while keeping the other end on the plate. 

To test that sidewalk have negligible effect on the dy- 
namics, we also performed experiments using a container 
with a shallow 1mm deep U-shaped groove with crossec- 
tional curvature which is at least four times that of the 
spheres in the dimer. The curvature focuses the dynam- 
ics to a vertical plane without dimer-sidewall collisions. 
Similar results were observed. However, higher quality 
position data is obtained with the container with side- 
walls, and used to report details of dynamics. 

We can characterize the different modes by the total 
mechanical energies of the vibrating dimer. As a simple 
measure of such energy we take the sum of maximum 
elevations of both spheres above the plate. Ignoring in- 
elasticity and horizontal translation yields the following 
estimates of energies per unit mass of drift, jump, and 
flutter modes E d = gT 2 /8, Ej = gT 2 /A, E f = gT 2 for 
A r 3> 1. Therefore the ratio of the energy for the drift, 
jump and flutter mode is approximately 1:2:8, which is 
consistent with observations. 

In the remainder of the paper we focus on the D-mode, 
since it exhibits the most non-trivial dynamics. The hor- 
izontal coordinate of the left tip of the dimer "resting" on 
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FIG. 3: (a) Horizontal position of the left tip of the dimer as 
function of time in the drift mode for A r = 3.9 (solid line) 
and 5.7 (dashed line), (b) The corresponding horizontal tip 
velocity. The thick solid line indicates the vertical position of 
the plate (/ = 25 Hz, F = 0.9). 



the plate is plotted as shown in Fig. [3(a) versus time for 
dimer with A r — 3.9 (solid line), and A r — 5.7 (dashed 
line) at / = 25 Hz. In Fig. |Sfa), the net motion of the 
A r = 3.9 dimer is directed towards its bouncing end. In 
the following, this direction is called forward or positive. 

We measured the horizontal drift speed of dimers with 
various A r between 2 and 10 at three different frequencies 
of vibration [see Fig. Ufa)]. The speed is determined by 
measuring the time taken by the dimer to move over a 
55 mm horizontal distance. For A r < 2.5, the D-mode 
becomes unstable. The drift speed is maximal for an 
aspect ratio of approximately 3.5, and then observed to 
decrease as A r is increased. It can be also noted that 
dimer moves in the backward direction at large A r > 5. 

We also varied the driving parameters to check their 
influence on the dimer dynamics. We found that the drift 
speed of the D-mode is inversely proportional to / [see 
Fig. 0(a)]. However, for fixed /, the drift speed is almost 
independent of T over a broad range. For example, for 
A r = 3.9 and / = 25 Hz, the measured speed was approx- 
imately constant over a range of T between 0.5 and 1.1. 
For smaller T < 0.5, the D-mode spontaneously collapses 
to rest, and for higher T > 1.1, the D-mode is unstable 
with respect to the transition to the higher energy J and 
F-modes, and their combinations. 

We checked that the drift motion is not due to any 
residual tilt in the plate by testing that the speed is the 
same when the dimer moves in either direction. Further- 
more, when the oscillated plate was deliberately tilted, 
uphill motion was observed depending on A r over a range 
of angles [see Fig. B[b)] . (These experiments where per- 
formed with the U-shaped bottom plate which has no 
sidewalk.) Thus the driving mechanism in the drift mode 
is robust enough to overcome gravity. 

Detailed examination of each tip reveals that the hor- 
izontal drift occurs due to the asymmetric slip of the 
low-energy left sphere "resting" on the plate (see Fig.|3J). 
Specifically, the horizontal velocity Ui of the left tip is 
plotted in Fig.dtb). For small A r , the positive slip which 
begins immediately after the right sphere collides with 
the plate, gradually slows down and is succeeded by a 
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FIG. 4: (a) The mean horizontal speed («) of the dimer scaled 
by gravity and frequency versus A T for three values of driving 
frequency. Positive direction points towards the bouncing end 
of the dimer. Symbols corresponds to experiments, and the 
line to the event-driven simulations at / = 25Hz (lines cor- 
responding to / = 50, 75 are virtually indistiguishable from 
this line and are not shown). Parameters of simulations are 
r = 0.9, e = 0.65, fi = 0.2, ^ s = 0.4, and p, = 0.1, fl s = 0.2. 
(b) Phase diagram of the direction of dimer translation in 
the D-mode as a function of plate inclination angle and A r . 
Robust uphill motion is observed for a broad range of param- 
eters. 



slower negative slip during the phase when the right end 
of the dimer approaches the plate, after which the process 
repeats. For large A r , the positive slip is much smaller 
or absent, and the drift is dominated by negative slip. 

Although it appears from Fig.^c) that the left sphere 
is simply rolling on the plate, in fact it performs a se- 
quence of small jumps immediately after the right sphere 
collides with the plate. We confirmed both the loss of 
contact after the collision during the small jumps as well 
as the fact that the left sphere is in contact with the plate 
for a large fraction of the cycle by measuring the electric 
resistance between the spheres and the plate. 

Our theoretical description of the dimer dynamics is 
based on a detailed analysis of frictional contacts between 
the dimer and the plate. Frictional inelastic oblique col- 
lision between asymmetric bodies has been a subject of 
many studies in contact mechanics where various re gim es 
of isolated contacts have been identified (see e.g. [ill]). 
Depending on the geometry and the incident velocity of 
the body one can distinguish three qualitatively different 
regimes: continuous slide, slip-stick, or slip-reversal (see 
also Ref. |5j for the case of thin rods). In addition to 
these regimes, the dimer in the D-mode experiences an 
even greater variety of contacts with the plate, includ- 
ing three types of double collisions (slide, slip-stick, and 
slip reversal) when two spheres touch the plate simul- 
taneously, and two regimes of rolling (with or without 
sliding). 

The dynamics of the dimer is governed by the Newton 
equations written for the center of mass (CM) velocity 
v = (u, 0, v) and the angular velocity u) 

mv = ^F c + mg, /u = ^r c xF c , (1) 



where to, I are the dimer mass and moment of iner- 
tia, F c = (F£,0,F£) is the contact force, and r c = 



(—X c ,0,—Z c ) is radius-vector from CM to a contact 
point (CP) [see also Fig.Q. The sum over c G {I, r} indi- 
cates that either left or right or both spherical ends can be 
in contact with the plate during collisions. The CP veloc- 
ity components in the plate frame of reference are given 
by U c = u — u)Z c , V c = v + ujX c — V p i , where V p i (t) is the 
plate velocity. In a sliding contact, F£ — —j^iF^, where fi 
is the friction coefficient and j = sgnf/ c , and in the stick- 
ing phases, F£ is found from the no-slip condition U c = 0. 
Transitions from stick to slip occur when \F£\ = ^ S F^ 
with the static friction coefficient fi s . Although the grav- 
itational force can be neglected during collisions, it must 
be taken into consideration during flight and rolling. 

For isolated collisions, we integrate Eqs. 0) with ini- 
tial conditions (u b ,v b ,ui b ) and the kinematic condition 
V c a = — eV b (where indices b, a correspond to before and 
after collision, and e is the kinematic restitution coeffi- 
cient), to obtain after-collision CM velocities (u a , v a , u> a ). 
Note that the dynamic and static friction coefficients for 
the rolling phase in general can differ significantly from 
those for the short collisions 0] and we use separate 
notations fi and fi s for the latter coefficients. A detailed 
description of all contact types is rather cumbersome and 
will be presented elsewhere (see also Supplementary 
material) . 

We combined the conditions for the occurrence of the 
various kinds of collisions and the forces acting on the 
dimer into an event-driven algorithm. We used this al- 
gorithm to elucidate the origin of the observed regimes 
of dimer motion. In D-mode, the left tip rolls on the 
plate with negative slippage as the right tip approaches 
the plate (see Fig. [SJ). At time to, the right tip col- 
lides with the plate with an incident contact velocity 
(U b < 0,V b < 0). In the first phase of collision, neg- 
ative slippage slows down due to the frictional force 
F? x = /1FJ. For large enough friction, at a certain 
time t' Q > to during collision, the negative slip stops 
(U r (t' ) = 0,V r (t' ) = V*). If formally V r s > ~eV r b , 
the dimer leaves the plate before the negative slip stops, 
and so the collision would be of a continuous slide type. 
Continuous negative slide cannot persist in the station- 
ary regime, because this would yield cumulative positive 
horizontal impulse transfer from the plate to the dimer 
over the period of vibration and hence an accelerated 
motion. For V" < —eV b , there can be two possible sce- 
narios for the second phase of the double collision. For 
large enough aspect ratios, fi s (A 2 + 7/5) > A, (where 
A = A r - 1,) the CP sticks at t > t' , i.e. U r (t > t' ) = 0. 
For jj s (A 2 +7/5) < A, the slip reverses direction, in which 
case it can be shown that is given by 



5jiA 2 - 5A + 7/t 
10A 2 - 5/L4 + 2 



(1 + e)V b 



10A 2 + 5/L4 + 2 
5fiA 2 + 5 A + 7/t 



U. 



For small A r , becomes positive after the double 
collision in agreement with experiment [cf. Fig.[SJb) and 
Fig. EHb)]. This "kick" during the two point collision is 
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FIG. 5: Vertical coordinates of left and right tips of the dimer 
(a) and the horizontal velocity of the left contact point for the 
short dimer A r = 3.9 (b) and long dimer A r = 6.1 (c). Differ- 
ent symbols and line styles denote different modes of motion 
and collisions between the dimer and the vibrating plate: sin- 
gle slide (x), single slip reversal (□), double collision (A), 
rolling (thick solid line), flight (thin solid line). Parameters 
of simulations are the same as in Figure 2] 



the main source of the positive drift momentum. Simul- 
taneously, the right tip acquires a finite vertical velocity 
which equals —eV^. The left tip also acquires a small 
positive vertical velocity <C , resulting in the dimer 
being lifted off the plate. Following the lift-off, the left 
end makes a sequence of collisions with the plate which 
result in a decaying sequence of small flights. This se- 
quence converges through inelastic collapse at time ti, 
after which the left tip of the dimer stays on the plate. 
Meanwhile the right end of the dimer is still in the air, 
and the rolling mode with slippage commences and con- 
tinues until the right tip finally collides with the plate, 
after which the process repeats. 

For large A r , the double collision is of the slip-stick 
type, so that Z7" = 0, and the mean slip velocity is nega- 
tive [Fig. Etc)]. While this picture qualitatively agrees 
with the experimental findings [compare Fig. 15(b), (c) 
with Fig. |3(b)], in order to match the mean drift veloc- 
ity quantitatively [see Fig. 0(a)], we had to choose p,, jx s 
for short-term collisions which were half of those for the 
rolling phase. While we do not have a good explanation 
for the smaller values, it is well known |l2j that the fric- 
tion properties depend strongly on the age of the contact 
and therefore can differ between short-time collisions and 
relatively long rolling contacts. As seen in Fig. 0(a), even 
with this fitting of friction parameters, there is a signifi- 
cant deviation between theory and observations for large 
aspect ratios. In particular, experiments do not show a 
discontinuous transition from forward to reverse drift at 
a certain A r , We believe that this discrepancy can be ex- 
plained by the sensitivity of the double collision to small 
variations in frictional properties of the plate at large A r . 



In conclusion, we investigated the dynamics of a dimer 
bouncing on an oscillated plate with the special focus on 
the first excited D-mode in which the dimer travels in 
one of the two directions selected via a symmetry break- 
ing. In typical fluctuation-driven transport systems, the 
direction of transport is selected by the local asymmetry 
of microscopic potential [||. In contrast, our bouncing 
dimer presents a novel phenomenon of ratchet motion in 
a system without an imposed asymmetric potential . 
Our theoretical analysis reveals that the D-mode is in fact 
a concatenation of multiple collisions which include indi- 
vidual, double, and rolling contacts. Horizontal trans- 
port is forced by the Coulomb friction during collisions 
and rolling phases. Because our preliminary experiments 
show that the drift mode survives collisions, one can an- 
ticipate that a collection of dimcrs, while initially cxhibt- 
ing random motion, will align mutually due to excluded 
volume interactions amplified by inelastic collisions. We 
expect that this alignment would give rise to a global 
convection similar to vortices observed with rods 0, E| • 
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